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We have used particle tracking methods to study the dynamics of individual balls comprising a 
granular flow in a small-angle two-dimensional funnel. We statistically analyze many ball trajectories 
to examine the mechanisms of shock propagation. In particular, we study the creation of, and 
interactions between, shock waves. We also investigate the role of granular temperature and draw 
parallels to traffic flow dynamics. 

PACS Numbers: 45.70.Mg 



I. INTRODUCTION 

Density waves are known to propagate in a variety of 
granular flow systems. For example, flows in funnels or 
hoppers can have slowly propagating density waves that 
may move up or down depending on sand properties and 
geometry 0. They also occur in long pipes M and closed 
hourglasses § as a result of the interaction with the in- 
terstitial fluid. In particular, they have been observed 
in a flow of monodisperse balls in both small-angle 
and wide-angle Q two-dimensional funnels. In the for- 
mer, one in fact observes kinematic shock waves which 
propagate upstream. 

The kinematic shock waves in a small-angle two- 
dimensional funnel have a complex phenomenology of 
their own and have been studied extensively in a pre- 
vious work (|] (hereafter referred to as HD). It was es- 
tablished that the flow can be divided into three different 
types depending on the funnel half-angle (3, each having 
a different characteristic shock speed U : 

Pipe flow ((3 < 0.1°): the shocks can be stationary 
or very slow (U < 10 cm/s) but always stop before 
reaching the inlet to the funnel. 

Intermittent flow (0.1° < (3 < 0.5°): strong- 
shocks propagate the full length of the funnel with 
50 < U < 150 cm/s. 

Dense flow (j3 > 0.5°): the flow is densely packed 
and the shocks are weak but fast (100 < U < 
300 cm/s). For (3 > 2.0°, shocks can no longer 
be observed by eye. 

In general, it was found that both the average local shock 
speed U(x) and shock frequency v{x) increased with the 
distance from the funnel outlet x. It also appeared that 
U(x) ~ w(x)/D where w(x) is the local funnel width and 
D is the funnel outlet width. Moreover, it was established 
that the shocks are mainly created at specific positions 
as a consequence of the monodispersity. Finally, different 
types of interactions between shocks were also observed. 

In the present work we simultaneously track all balls 
in the field of view of a camera to study the flow on 



smaller length and time scales. This enables us to observe 
shock creation and interaction more closely. We can then 
compute the granular temperature and also traffic flow 
curves. The paper is structured as follows. In Sec. || we 
describe the experimental setup and the particle tracking 
method. In Sec. Ill we study the properties of individual 
ball trajectories. In Sec. iTVl we reduce the ball flow to 



we discuss how 
we describe their 



various one-dimensional fields. In Sec 
shock waves are crea ted a nd in Sec. [VI 
interactions. In Sec. VII we dis cuss the role of granular 
temperature in the flow. In Sec. VIII we present the time 
averaged flow properties. In Sec. IX we compare the ball 
flow with traffic flow. Finally, we summarize our results 
in Sec. g 
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FIG. 1. (a) Schematic of the experiment, (b) Parameters 
of the funnel geometry. 



II. DATA ACQUISITION 



A. Experimental setup 



The experimental setup is described in detail in HD 
and will only be briefly reviewed here. A top view is 
shown in Fig. 0(a). The granular matter consists of 50000 
3.18 mm diameter brass balls. They roll in one layer 
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between 3.45 mm high aluminum walls (A) on a coated 
Lexan plane (B). The aluminum walls have 200 cm long 
straight sections which then open smoothly at the top 
forming a reservoir area (C). The straight sections can 
be moved to vary the outlet width D and the half-angle 
j3 of the funnel (see Fig. 0(b)). The Lexan plane is tilted 
an angle 8 so the balls flow through the funnel into a 
collection container (E). Another plate (not shown) was 
placed on top of the aluminum walls to keep the flow in 
a single layer. Unless stated otherwise, the outlet width 
and inclination angle were kept fixed at D = 10 mm and 
8 = 4.1°, respectively. The funnel half- angle was varied 
in the range /3 = 0° - 3°. 

Elements of the data acquisition system are also shown 
in Fig. |l|(a). A light box (F) was placed below the lower 
120 cm of the funnel and a video camera (G) was placed 
above it. Its analog output signal (H) was sent to a frame 
grabber in a PC (not shown). The camera was mounted 
on a stiff bar (I) and could be moved along its length. 
The bar was supported by a stabilized stand (J). 



B. Video system 



The video camera is a Pulnix TM-6701AN, 8-bit grey 
tone, non-interlaced, analog CCD camera. It is capable 
of filming in four modes of which we used two: 640 x 200 
pixels at 130 frames/second (fr/s) and 640 x 100 pixels 
at 221 fr/s. The camera shutter may be adjusted from 
no shutter to 1/32000 s. Shutter speeds of 0.5-1 ms are 
mostly used due to the available light intensity. The ana- 
log signal from the camera is not a standard video signal 
but has a pixel rate of 25.5 MHz. This signal is trans- 
mitted to a Matrix Vision PCimage SGVS frame grab- 
ber card. This PCI-bus PC-card has almost no internal 
memory and data is transmitted by direct memory ac- 
cess to the main memory of the host computer. Up to 
32 MBytes may be grabbed per measurement. Thus, for 
(3 < 0.1°, up to 2000 consecutive frames may be grabbed. 
For the largest angles only ~ 400 consecutive frames can 
be stored. Consequently, recorded sequences typically 
span 2-9 s, corresponding to 1-15 shocks depending on j3 
and x. 

The camera is equipped with a 16 mm C-mount lens 
with a field of view of ~ 20°. The camera is oriented such 
that the 640 pixel direction is parallel with the center 
line of the funnel. The view length is used to denote 
the length of the funnel visible in these 640 pixels. The 
lens was set at a height of 106 cm giving a 36.9 cm view 
length. The camera was placed so as to cover one of the 
following ranges: x — —0.2-36.7 cm, x — 36.8-73.7 cm, 
ori = 73.8-110.7 cm. 
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FIG. 2. Sequence of video frames showing a propagating 
shock (/? = 0.15°, D — 10 mm, 6 = 4.1°). Only every fourth 
frame is shown. 

The aluminum walls are painted black and extended 
with black plates to ensure that the only non-black part 
of the camera's field of view is the funnel lighted from 
below. The balls therefore appear as black "dots" in a 
white funnel as shown in the frame sequence in Fig. 0. 



C. Ball tracking method 

At 221 fr/s, a ball speed of 1 ball diameter per frame 
corresponds to a velocity of ~ 70 cm/s. Dense flows must 
be slower to be tracked. In low density flows, tracking is 
possible at significantly higher speeds due to the larger 
ball separations. (Computer based particle tracking is 
the only feasible method for data sets containing up to a 
million ball images.) 

Previously, several experiments with fast, low den- 
sity flows have been recorded on high speed film (~ 
1400 fr/s), with ball positions then obtained by eye 0. 
This method is useful for following relatively small groups 
of particles over relatively short periods of time. Warr 
et al. H have used expensive digital high speed video 
equipment to follow a few (27-90) fast particles. Particle 
positions are subsequently determined by computer using 
the Hough transformation (an edge detection method ) 
for determining ball perimeters. This method requires 
many pixels 150) for each ball. By comparison, we 
use ~ 30 pixels/ball. 

We now sketch the algorithm used in the present exper- 
iment. (The full details are described in Appendices |A] 
and ||.) Each frame is subtracted from the image of 
the empty funnel. Thus, the balls stand out as a white 
patch on a black background. A threshold is then de- 
termined based on local grey scale distributions and ev- 
erything above the threshold is then defined as part of 
a ball. Each ball center position and height (i.e., grey 
scale intensity) is estimated, and a number of pixels are 
assigned to the ball image (typically 20-35). These data 
are used as input to a fitting routine. The pixel values 
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(with coordinates (x,y)) are fitted to a two-dimensional 
peak function and the center of the fitted peak position 
determines the ball center (x J c (tk),yi(tk)) for the jth ball 
in the fcth frame. The final result is illustrated in Fig. [|. 
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FIG. 3. Comparison of (a) video frames and (b) cor- 
responding artificial frames using ball tracking (/3 = 0.3°, 
D = 10 mm, and 6 — 4.1°). In (b) each ball is identified by a 
shade of grey. 

Each frame has now been converted into a table of ball 
positions. Based on the correlations between the first two 
frames and their densities, an initial guess is made of the 
one-dimensional coarse-grained velocity field v{x,t\) of 
the first frame. The jth ball in the first frame is now 
assigned a velocity = ( v [ x i( t i)^ t i]^)- I n 

the interest of clarity, we will henceforth drop the in- 
dices j and k. Based on these velocities, predictions of 
the possible positions in the next frame are determined 
and compared with the actual ball positions. Positions 
are matched "one-to-one" making the best matches first. 
Balls in the old frame that cannot be matched accept- 
ably are defined as lost. Matching is first attempted at 
the position predicted for constant (v x (t), V y (t)), then 
at the old position (v x (t),v y (t)) = (0,0), and finally 
at a midway point. Unmatched new balls are assigned 
(v x (t),v y (t) — (v(x, t),0) at their respective positions. 
(In each frame v(x,t) is re-computed from averages of 
those balls which have known velocities.) By repeated 
application of this method, the full trajectory of each 



ball can be found. 

The error of (x c ,y c ) is < 0.15 mm, but in a few bad 
cases it may be as much as ~ 0.5 mm. Less than 0.05% 
of the balls are normally lost in each frame. Most balls 
can be followed from the time they enter the view length 
of the camera until they leave it. Bad matches are very 
rare. (They only seem to happen in groups and are easily 
recognized.) 
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FIG. 4. Histograms of measured ball center positions, (a) 
f3 = 0.1°. One only observes weak structures, (b) /3 = 0.5°. 
There are now 4 — 10 triangular packing rows, (c) /3 = 1.5°. 
There are 4—17 triangular packing rows. (The vertical white 
stripes separate independent data sets. The weak horizontal 
stripes in (a) and (b) with a separation of ~ 0.6 mm are 
caused by interference with camera pixel rows. See Appendix 
B.) 

The final determination of the velocity (v x (t),v y (t)) 
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and acceleration (a x (t),a y (t)) of each ball is based on 
a least squares fit of the positions to a parabola using a 
number of consecutive points. This number represents a 
tradeoff between precision and time resolution. For ve- 
locities, 5 points are generally used, resulting in a time 
resolution of ~ 23 ms. For accelerations, 5 or 7 points 
are used depending on the application. The latter corre- 
sponds to a time resolution of ~ 32 ms. 



III. FLOW IN LAGRANGIAN COORDINATES 



A. Single ball trajectories 



Since each measurement typically involves 0.3-1.0 mil- 
lion balls, it is possible to measure high resolution distri- 
butions of the ball center positions (x c (t),y c (t)) (averaged 
over all frames). Some examples are shown in Fig. [jj For 
intermittent flow (Fig. ||(a)) we observe weak periodic 
patterns across the funnel. For dense flows (Figs. f|(b) 
and (c)), these patterns are much more pronounced. (No 
such patterns are observed for pipe flow.) 

As discussed in HD, the monodispersity of the balls 
allows triangular close packing at packing sites x = Xi 
given by 



It + V3r(i - 1) - D 
2tan/3 



(1) 



where i is an integer and r is the ball radius. This pre- 
dicts that such packing sites should recur at intervals of 
y/3r and v / 3r/2tan/3 in the funnel axis and transverse 
directions, respectively. This is what is indeed observed 
as can be seen in Figs. |(b) and (c) to within - 3%. Im- 
mediately downstream of a packing site the lattice must 
rearrange itself. It appears that this occurs over a rel- 
atively short distance. For example, in Fig. |^(b) where 
X7 = 56 cm and xs = 72 cm, the rearrangement oc- 
curs at x — 70 cm. As one might intuitively suspect, 
shock waves are more readily produced where collisions 
are most likely. We believe this is the reason for the 
nearly exclusive creation of shock waves at the packing 
sites Xi as discussed in HD. 

Fig. |5|(a) shows four representative examples of individ- 
ual ball traj ectories obtained using the method described 
in Sec. |ll C| . Figs. ||(b) and (c) show v x (t) for the same 
four balls. Note that since the balls almost always move 
downstream, we always plot — v x (t) for convenience. It 
is apparent from these figures that an individual ball ex- 
periences long periods of constant (negative) acceleration 
interrupted by abrupt drops in speed to \v x \ ~ 1 — 5 cm/s 
for (3 > 0.2° and ~ 5-15 cm/s for (3 < 0.1°. These drops 
occur when the ball encounters a shock. This qualitative 
behavior is even more apparent in the acceleration a x (t) . 
In Fig. p] we show the velocity v x (t) (a) and acceleration 
a x (t) (b) for another trajectory. 




FIG. 5. (a) The trajectories of four balls for different 
funnel angles. The a; c (t)-coordinate is plotted for each ball, 
(b) The velocity v x (t) of the two lower trajectories in (a) 
(small funnel angles), (c) The velocity v x (t) of the two upper 
trajectories in (a) (large funnel angles). 



B. Collisions and the coefficient of restitution 

By observing collisions between two balls, we can mea- 
sure (with limited resolution) the coefficient of restitution 
e of a ball, defined as the ratio of the relative velocities 
before and after a collision. This can only be done for 
flows with low densities and no shock waves in the field 
of view so that individual collisions can be identified. In 
two data sets at f3 = 0° that adhere to these require- 
ments, we have measured a total of 192 collisions and 
found that e = 0.78 ± 0.02. The average relative velocity 
before these collisions was v re i — 7 ± 3 cm/s, compared 
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with the actual ball speed = 112 ± 6 cm/s. Assuming 
that the balls roll without slipping, the tangential veloc- 
ity difference immediately preceding each collision should 
be twice the ball speed and therefore much higher than 
v re i. The measured value of e using this method agrees 
with the value e > 0.74 found by dropping a ball on a 
hardened steel block. 

From measurements of fast rolling collisions (v x ~ 40- 
80 cm/s) with a steel barrier, we estimate e ~ 0.5- 
0.6. The additional energy loss here may come from the 
rolling ball sliding on the contact surfaces at the moment 
of collision. This may be related to the situation where 
a fast ball encounters the slowly moving packed region of 
a shock wave. 



C. Shock profiles 

We showed in Figs. ||(b,c) and ||(a) the velocity v x {t) 
of single ball trajectories. Their apparent sawtooth be- 
havior is, in fact, a general feature observed in many 
trajectories. From this simple structure we can obtain 
some interesting information concerning the dissipation 
of the flow between shock waves. 
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FIG. 6. Movement of a single ball at j3 = 0.8° between 
x = 73 cm and x — 37 cm showing the (a) velocity v x (t) 
and (b) acceleration a x (t). The peaks in (b) indicate that the 
ball has encountered and is nearly stopped by a shock. The 
dashed line marks zero acceleration. 

First, we determine the inter-shock acceleration a, s , 
i.e., the acceleration between shock waves. Using the 
acceleration a x (t) shown in Fig. ^(b) and employing a 
threshold criterion, we can identify shock waves and the 
time t s when a given ball encounters a shock wave. Lin- 
ear fits yield typical values of a is = —40 ± 10 cm/s 2 . 
This value should be compared with the expected effec- 
tive gravitational acceleration (for spheres rolling with- 
out slipping) g r = If? sin = — 50 cm/s 2 for 9 = 4.1°. It 
is not surprising that \ai s \ < \g r \ since energy is probably 
lost in collisions and friction processes. 




FIG. 7. Average velocity profiles of shocks in Lagrangian 
coordinates. Specifically, we plot the average of —v x (t — t 3 ) 
of individual balls, where t s is a pronounced local maximum 
in ball acceleration, (a) (3 = 0.15°. Average of one shock, (b) 
P — 2.0°. Average of ~ 10 shocks, (c) single shock profile 
for (3 = 0.2° and different inclination angles 6. The dashed 
lines in (a) and (b) indicate the standard deviation of the 
measurement. 



Better statistics can be obtained by averaging over 
many trajectories of balls encountering the same shock 
wave. By using a film sequence such as that in Fig. ^, we 
choose a time interval such that only one shock is visible 
in the view length. Single ball trajectories with values 
of t s within that interval are then selected (typically a 
few hundred). If we now average the trajectories for all 
these balls by setting t s = for each ball, we obtain 
average profiles such as those shown in Figs. |?](a) and 
(b) for intermittent and dense flow, respectively. These 
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preserve the sawtooth structure discussed earlier, rein- 
forcing our picture of the behavior of single balls passing 
through shock waves. In particular, the inter-shock ac- 
celeration calculated for points outside the shock region 
for 8 = 4.1° is a is = —38 ± 3 cm/s 2 . It was found to be 
independent of both velocity (i.e., the same just preced- 
ing or following a shock wave) and (3 when (3 > 0.1° (i.e., 
for intermittent and dense flows; there are exceptions in 
pipe flow when \v x \ > 80 cm/s). It was also checked that 
ai s was the same in different regions of the funnel, and 
hence independent of x as it should be. 




w/D 

FIG. 8. (a) The median of the distribution of v m in versus 
x. (b) The same data replotted against the rescaled variable 
w{x)/D. 

We have also checked the dependence of a,i S on 8. The 
velocity profiles for three different values of 8 are shown 
in Fig. |t](c). Ignoring for now differences in the velocity 
jumps across the shock, one sees that ai S increases with 
8. For 8 = 3.0°, we find a ls = -27 ± 3 cm/s 2 and for 
8 = 5.0°, ai S = —45 ± 3 cm/s 2 . (These values were also 
independent of velocity and (3.) It appears therefore that 
&is/9r ~ 0.75 in all three cases, implying that can be 
regarded as an effective gravitational acceleration. 

Experiments involving only a single ball for 8 = 4.1° 
yielded a measured value g r — — 45±5 cm/s 2 . The reason 
for the difference between this and the theoretical value 
g r = —50 cm/s 2 is partially due to a local deviation of 
8 of ~ 0.2° — 0.3° from a small bending of the plane, 
and possibly also to non-ideal rolling. Consequently, we 
believe that the difference in ai S and the measured value 



of g r is significant, i.e., due to the interaction between 
balls in intermittent and dense flows. 

Just before and after a shock, a ball reaches its max- 
imum and minimum velocity v max and v m in, respec- 
tively (see, for example, Fig. [fj]). In particular, knowl- 
edge of v m i n is necessary to understand the creation of 
shock waves discussed in Sec. |VT| , In the same manner 
as discussed earlier, we find that the average value of 
Vmin f° r au balls passing through a given shock wave is 
bmm| ~ 10 — 20 cm/s, 4—10 cm/s, and 1 — 5 cm/s, in 
pipe, intermittent, and dense flow, respectively (see for 
example Figs. [|, and [?]). 

We also checked the dependence of v m in on [3 and x. 
(Due to the presence of shock waves that are created 
within the view length and a low number of shock waves, 
v m in{x) is obtained as the median of the distribution of 
Vmin, not the mean.) This is shown in Fig. 0(a) for 0.4° < 
(3 < 3.0° and D = 7, 10, and 14 mm. (There were 
too few shock waves for smaller (3 to obtain meaningful 
statistics.) In HD it was found that the local shock speed 
U(x) depended linearly on the local funnel width w(x). 
If we plot Vmin against w(x)/D, we obtain the plot shown 
in Fig. ||(b). The apparent collapse of the data indicates 
that there may indeed be some simple scaling properties 
inherent in the system. It will be shown in a future work 
how the behavior of ai S and v m in can be used to model 



shock wave statistics |T(| . 



IV. FLOW IN EULERIAN COORDINATES 

A. Granular fields 

From the full set of ball trajectories in the flow, we 
can construct continuous one-dimensional Eulcrian fields, 
such as the density p{x, t) and velocity v(x, t). (The justi- 
fication for considering a one-dimensional flow is given in 
Appendix |c"[) Of course, to obtain meaningful results, it 
is necessary to coarse-grain. In order to obtain a spatial 
resolution better than one ball diameter, this was done 
by drawing a primitive unit cell (for close-packed cir- 
cles) around the ball center positions (x c (t),y c (tj), with 
two sides parallel to x (i.e., the center line of the fun- 
nel). The funnel was then binned into strips of width 
~ 0.1 mm and the occupied area was calculated for each 
strip. These strips can then be summed to achieve any 
desired resolution, although typically we used a resolu- 
tion of ~ 1 mm. We then divided by the area of the 
strip to obtain the relative density p(x, t). Thus, for close 
packed balls, p(x,t) = 1. (The relative density here dif- 
fers from that in HD which was measured using relative 
light intensities.) The velocity field v(x,t) is then com- 
puted as the weighted average 
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FIG. 9. The different fields for (3 = 0.8 
granular temperature T(x,i). 
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a) relative density p(x,t), (b) velocity v(x,t), (c) positive acceleration a(x,t), (d) 



where the sum is over all balls whose unit cells con- 
tribute to a strip and p 3 is the corresponding relative 
density. The acceleration field a(x, t) is obtained in a 
similar manner. 

The granular temperature per ball is defined as 
T(x, t) = ((v y -Vy) 2 ) + ((v x -v x ) 2 } where v x = v(x, t) and 
v y = 0, and where the brackets indicate an average as in 
Eq. ^. It can be written as the two components T(x, t) — 
T x (x,t) +T y (x,t), where T x (x, t) = ([v x - v(x, t)] 2 ) and 
T y (x,t) — (v^), in order to see if the granular tempera- 
ture is isotropic. 



ing approximately to the packing sites xs = 45 cm and 
X9 = 55 cm (see Eq. [l]). 

The granular temperature T(x,t) shown in Fig. ^(d) 
clearly peaks in the shock regions as it did for a{x, t), 
although it is not as evenly distributed along the shocks. 
The temperature appears to cool down over relatively 
short time and length scales. There is a background level 
of ~ 2 cm 2 /s 2 with a weak packing site periodicity of 
~ 10 cm. The granular temperature will be discussed 



more thoroughly in Sec. VII 



B. Discussion of the fields 

An example of the relative density is shown in Fig. ||(a) 
for dense flow (/? = 0.8°). It reveals five faint shocks. 
However, in dense flows, density fluctuations on the scale 
of a ball diameter tend to be comparable with or even 
dominate the fluctuations associated with the passing of 
shock waves. (For example, the whiter tracks moving 
downstream are defects, which are often vacancies. In 
dense flows they can survive the passing of shock waves 
since shocks do not seriously rearrange the local packing.) 

The corresponding velocity and acceleration fields are 
shown in Figs, ^(b) and (c), respectively. In the velocity 
field, the same shock waves are much more clearly visible 
since the greatest contrast occurs when fast balls enter 
a shock wave and lose most of their speed. (Note that 
again, since v(x, t) is always negative, we plot — v(x, t) 
here and in all subsequent plots.) 

By adjusting the resolution in the positive accelera- 
tion field a(x,t) shown in Fig. ^(c), we can improve the 
contrast in v(x,t) in order to reveal weaker details of 
the flow. For example, one can now see the creation of 
(temporarily) stationary shocks at locations correspond- 



C. Flow dynamics for fixed x 

The grey scale plots of the granular fields in Fig. || 
do not convey the magnitudes of the fields very well. 
Therefore we show v(x,t) and p(x,t) for four fixed x in 
Fig. [h] and study how these quantities are affected by 
propagating shock waves. 

In Fig. [H](a), the velocity v(x,t) is shown at x — 44 cm 
(midway between Xi = 24 cm and \5 — 64 cm), and at 
x = \b — 64 cm for an intermittent flow ((3 = 0.2°). 
(These positions were chosen for reasons which will be- 
come clear later.) One can clearly see propagating shock 
waves and smoothly increasing flow speeds between the 
shocks in the range 3-45 cm/s. From the displacements 
of the shock fronts, we estimate that U ~ 120 cm/s. At 
x = 64 cm, one can see two shocks labelled Al and A2. 
From fields such as those in Fig. ||, we know that these 
are newly created shocks. At x = 44 cm, disturbances are 
observed at A3 and A4 as the velocity increases between 
shocks. These are caused by balls that have passed the 
newly created shocks at Al and A2. The corresponding 
relative densities are shown in Figs. [l0](b,c). They reveal 
that the shock fronts are relatively sharp, but otherwise 
details of the flow are generally not as clear as in v(x, t). 
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FIG. 10. Time series of —v(x,t) and p(x,t) at fixed x showing propagating shocks for two different /3. (a) velocity for 
f3 = 0.2° at x = 44 cm (between X4 and xs) an d x = \5 — 64 cm. Al and A2 indicate newly created shocks, and A3 and A4 
show, respectively, the effect of these shocks 20 cm downstream, (b) density corresponding to (a) at x = 44 cm. (c) density 
corresponding to (a) at x = 64 cm. (d) velocity for f3 = 2.0° at x = 4.5 cm (between packing sites) and x — 26.5 cm (at a 
packing site), (e) density corresponding to (d) at x — 4.5 cm. (f) density corresponding to (d) at x = 26.5 cm. 
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FIG. 11. Acceleration field a(x,i) showing weak shock waves for large values of f3: (a) j3 = 1.0°, (b) /3 = 2.0°, (c) = 3.0°. 
The grey scale for the magnitude is shown at the top of each plot. 



In Fig. pJj|(d), we show v(x, t) at x = 4.5 cm (between 
X4 = 2.4 cm and \5 = 6-5 cm), and x = xio = 26.5 cm 
for a dense flow {(3 — 2.0°). Many shocks with speeds 
U ~ 350 cm/s are now visible. At x = 4.5 cm they 
are all newly created. These shocks can still be observed 
at x = 26.5 cm along with many other shocks created 
upstream of them. The corresponding relative densities 
are shown in Figs. |l(](e) and (f), respectively These are 



of limited value when f3 > 1.5° Near the outlet (e) there 
are still fluctuations related to shock waves, but upstream 
(f) the density is constant and nearly unity despite the 
presence of many shock waves. 
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D. Shock behavior for /J > 1.0° 



As (3 is increased, the shock waves are faster, weaker, 
and more frequent as discussed in HD. For f3 > 1.5°, one 
can no longer study shock waves from directly measured 
density fluctuations as in HD and in this situation the 
ball tracking method is particularly useful. 

In Figs. |ll|(a-c) we show a(x, t) near the outlet for 
three relatively large values of f3. One can see that as 
(3 increases, the frequency and speed of the shocks in- 
crease while their strength decreases. In particular, in 
Fig. 11(c) it also appears that the shock waves die off 
noticeably with increasing x. 
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FIG. 12. The spatial dependence of (a) the velocity 
—v(x,t) and (b) the corresponding relative densities p(x,t) 
for fixed t. In these data, one is observing a decaying propa- 
gating shock. 
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FIG. 13. (a) Estimated shock speed U(x) for 

1.0° < (3 < 3.0° and D = 7, 10, and 14 mm. (b) The same 
data replotted against the rescaled variable w(x)/D. 



As already discussed in Sec. Ill C , it was observed in 
HD that for 0.1° < (3 < 1.0° the local average shock speed 
seemed to depend linearly on the local funnel width. We 
can now test this hypothesis for larger (3. Unfortunately, 
with our current data constraints, we have only 10-20 
shocks so the statistics are not as good as could be de- 
sired. Fig. |l3| shows U(x) based on a simple ridge de- 
tection method applied to a(x, t). Fig. [l3](a) shows U(x) 
for different values of (3 and D while (b) shows the same 
data versus w(x)/D. The data collapse in Fig. |l3|(b) 
is not complete, although it is consistent with the same 
analysis in HD. 



Once again it it useful to look at the magnitudes of 
v(x,t) and p(x,t), but now for fixed t. In Fig. ^ we 
show v(x,t) and p(x,t) for four equally spaced consecu- 
tive times. In Fig. |l2j(a) one can see a shock propagating 
upstream with a speed U ~ 250 cm/s. As already noted, 
it dies off as it moves upstream, but now one can clearly 
see that the shock region also gets broader. Fig. |l2|(b) 
indicates why the density p(x, t) cannot be used to de- 
tect shock waves for large j3. Except for the lowest 10 cm 
of the funnel, the density fluctuations are entirely dom- 
inated by the packing sites (which in this case have a 
periodicity of 3.2 cm). 



V. SHOCK CREATION 

Having established some understanding of how shock 
waves propagate, we would now like to examine under 
what circumstances they are created. Since for f3 > 0.1° 
shock waves (once created) propagate all the way to the 
reservoir, the mechanisms of shock creation will strongly 
effect the flow properties. In HD it was found that for 
intermittent flows new shock waves were almost exclu- 
sively created at the packing sites Xi- Fig- @ illustrates 
how a moderately dense flow is forced to reorganize at 







the packing sites, thus possibly increasing the likelihood 
of congestion there. 
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FIG. 14. Examples of starting shocks and their respec- 
tive v m in(x) and Vmax(x) for different funnel angles: (a,b) 
/3 = 0.2°, (c,d) (3 = 0.6°, (e,f) = 1.5°. 



In Sec. HI C it was discussed how v m i n varies in the 
different flow types. Basically v m i n is a measure of how 
effectively the shock has absorbed momentum and energy 
from the flow upstream of the shock. Generally, we have 
found that v min for newly created shocks is often signifi- 



cantly higher than the typical values of v m i n (x) shown in 
Fig. 10(a) where the shocks at Al and A2 are known to 
be newly created. This means that the new shocks do not 
dissipate energy as efficiently as the older shocks around 
them. Three examples of newly created shocks are shown 
in Figs. [u](a,c,e), and the averages of v m i n and v max of 
the balls encountering these shocks are shown in (b,d,f), 
respectively. One can observe how v m i n starts relatively 
high, but then over a distance of 5-15 cm, it drops to the 



typical values found in Sec. IIIC 
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FIG. 15. Time series of (a) —v(x,t) and (b) p(x,t) at 
x = 30 cm showing the creation of a shock in pipe flow 
fj3 = 0.075°). A slow, dense group of balls passes down- 
stream (Bl) and forms a shock at x — 26 cm (not shown) 
which starts propagating upstream (B2). 

For each flow type we can make more specific obser- 
vations about how shocks are typically created. In pipe 
flow new shocks typically start as a group of balls mov- 
ing downstream which have a somewhat smaller veloc- 
ity and higher density than those around them. This 
early stage of a shock at some point becomes stronger, 
stops, and then starts moving upstream as a stable shock 
(shocks moving downstream are generally not very sta- 
ble). An example of this process is shown in Fig. [l5| 
where a disturbance passes the observation point at Bl, 
turns around downstream, and then starts propagating 
upstream, passing the observation point again at B2. 
Mass conservation across the shock discontinuity requires 
that 



U = (p u Vu - PdVd)/{Pu - Pd) 



(3) 



where the subscripts u and d signify upstream and down- 
stream of the shock, respectively. In this example, 
p u ~ 0.3w, pd ~ 0.6w — + 0.8u>, and v u ~ —95 cm/s. 
Note that as Vd increases from —70 cm/s to —20 cm/s, 
and pd increases slightly, U changes sign. 

In intermittent flow the shock may start in the manner 
described above for pipe flow, but more commonly it is 
created at a packing site without the early downstream 
phase discussed above. An example of such a shock is 
shown in Fig. |l4|(a). 
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FIG. 16. Positive acceleration field a(x,t) showing examples of shock interactions, (a) Shock start repulsion (PO and PI) 
(/3 = 0.2°). (b) Two temporary pre-shocks (P2 and P3) (/3 = 0.4°). (c) Stationary shocks joining at packing sites (P4 and P5), 
and a near shock repulsion (P6), pre-shock (P7), and micro-shock (P8) (/3 = 0.8°). (d) Near shock repulsion (P9) (/3 = 1.5°). 



In Fig. 0(b) we see how the growing difference be- 
tween v m i n and v m ax shows how shocks improve their 
ability to absorb momentum and energy during the first 
10 cm of propagation. Other examples of shock cre- 
ation in intermittent flows have already been shown in 
Figs. |(a) and 0(a). 

In dense flow shocks are often created by the joining 
of a few relatively stationary weak "pre-shocks" (which 
will be discussed in the next section) or a region with 
seemingly "unstructured" disturbances in the flow at a 
packing site (see, e.g., Figs. || and 0). The nature of 
these shock creation processes is not known. Fig. |l4](c) 
shows the creation of a shock wave 38 cm upstream of 
the outlet. The |w m i„| shown in Fig. 0(d) initially has 
a value somewhat higher than what is typical for that 
funnel position (see Fig. ||(a)). In Fig. 0(e) the cre- 
ation of a shock near the outlet is shown for large (3. 
Fig. ^(f) shows how its |u m in| a ls° quickly drops after 
the creation. Near the outlet (and one packing site up) 
all shocks have just been created. Therefore, the v m m 
curve in Fig. 0(f) partially corresponds to the first part 
of the curves in Fig. ^. This also explains why the values 
of Vmin in Fig. 0(d) at x — 4.5 cm are relatively high. 



When the balls leaving this shock reach the next shock, 
they cause it to temporarily slow down (PI), in effect 
repelling it. (A similar patch of slower balls was observed 
in Fig. 0(a) at A3 and A4 following the creation of new 
shocks at Al and A2.) 

In Fig. 0(b) examples of pre-shocks (P2 and P3) sep- 
arated from the main shock by less than 0.1 s are shown. 
At P3 the pre-shock seems to gain amplitude and swal- 
lows the original main shock. Fig. 0(c) shows weak sta- 
tionary shocks (P4 and P5) waiting to be swallowed by 
arriving shocks at or just upstream of packing sites. Ex- 
amples are shown of a near shock repulsion at P6 and of 
a weak pre-shock at P7. Some micro-shocks also seem to 
propagate upwards with speeds similar to the surround- 
ing shocks at P8. It is not clear what these represent, but 
since the balls are densely packed at (3 = 0.8°, they may 
be "chain collisions" propagating in single rows of balls. 
The data set in Fig. |l^(c) partially overlap the data pre- 
sented in Fig. ^|. In Fig. 16(d) a near shock repulsion is 
seen at P9 and several weak waiting shocks are observed. 
These waiting shocks seem related to the shock creation 
process in dense flows discussed in Sec. |y|. There is no 
indication that shock interactions disappear at larger (3 
(compare with Fig. O). 



VI. SHOCK INTERACTIONS 

Shock wave interactions have been discussed in HD. 
Since these were studied using directly measured density 
fluctuations, they were hard to resolve and the nature 
of the interactions was somewhat speculative. The supe- 
rior spatial and temporal resolution of the ball tracking 
method now gives us a clearer picture of the interactions 
and even reveals new interactions not observed earlier. 

An example of a weak shock start repulsion is shown 
in Fig. 0(a). First, a new shock is created at point (P0). 



VII. SHOCK TEMPERATURE 



When a ball encounters a shock, it quickly loses most 
of its velocity (and therefore most of its momentum and 
energy) in collisions with other balls and the funnel walls. 
In the dense region of the shock where the mean distance 
between balls is small, this results in relatively high colli- 
sion rates when they encounter a fast incoming ball. The 
high density also forces balls to come into contact with 
the funnel walls so that the total momentum/energy of 
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a dense group of balls is also reduced by friction. This 
phenomenon is reminiscent of the clustering |p"l|| and in- 
elastic collapse |l2j observed in simulations of granular 
materials. 

This process can also be described as follows. Up- 
stream of the shock there is a steady flow of relatively 
high energy balls. When a group of balls encounters 
the shock, collisions "randomize" their translational en- 
ergy, increasing the granular temperature. However, it 
is then reduced by each collision and by the friction pro- 
cesses mentioned above. Some distance downstream of 
the shock the relative motion of neighboring balls appears 
to be negligible. Thus, since we cannot measure collision 
rates in a shock, the rate at which the granular temper- 
ature decays in a group of balls following the passage of 
a shock front seems the best way of studying the energy 
loss processes in the shock. Consequently, we define tt 
as the characteristic decay time of T(x, t) following the 
passage of a shock front in Lagrangian coordinates. 
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FIG. 17. Comparison of (a) T x (x,t) and (b) T y (x,t) show- 
ing creation and propagation of shocks for /3 = 0.6°. 

First, we must consider the usefulness of T(x, t) in the 
shock region. It may be argued that since v{x, t) is chang- 
ing rapidly in the shock region, a calculation of T x (x, t) 
(and thus T(x,t)) is of limited value here. The coarse- 
grain resolution on which T x (x,t) is based is < 1 ball 
diameter, while the shock region is typically > 3 ball 
diameters. Thus the resolution of T x (x,t) should suf- 
fice. In Fig. [l7] we compare T x (x, t) and T y (x, t) and find 
that in general they have the same magnitude, except 
for some small packing site related oscillations. (There is 
obviously less room for transverse movement of individ- 
ual balls near packing sites.) Consequently, we will only 
consider T(x, t) rather than the individual components. 

Transforming T(x, t) to Lagrangian coordinates to 
study the magnitude of tt would be difficult. Instead 



the following observation can be made. The character- 
istic decay length Ay of T(x, t) downstream of the shock 
can be estimated from graphs of T(x, t) for fixed t. Fol- 
lowing the passage of a shock front by a group of balls, 
the shock will move upstream with a speed U while the 
balls move downstream with a speed ~ |v m m|- Conse- 
quently, we find that tt ~ Xt/{U + \v m in\)- 
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FIG. 18. Spatial dependence of T(x,t) at the position of 
a propagating shock at subsequent times. In (a) the shock is 
created at x = 40 cm (xs = 41 cm) just before to- /.From 
to + 0.1 s to to + 0.35 s the shock propagates from x = 44 cm 
to x = 70 cm with increasing speed. In (b) the passage of the 
shock at the next packing site (xe, = 67 cm) is shown with 
better time resolution. 

The spatial dependence of T(x, t) around a propagat- 
ing shock is shown in Fig. [l^ for a series of fixed times. 
In Fig. |l8|(a) a shock wave is created at X5 = 41 cm and 
propagates from x = 40 cm to x — 70 cm in 0.35 s. The 
propagation past xe — 67 cm is shown with higher tem- 
poral resolution in Fig. |l8|(b). Due to the dense packing 
at and just below a packing site (62 cm < x < 67 cm in 
Fig. |l|(b)), there is little room for relative motion be- 
tween the balls. When a shock passes this region, T(x, t) 
remains small but has peaks both upstream and down- 
stream. The tendency of T{x 1 1) to remain small in ar- 
eas just below packing sites is generally observed (see 
Figs. |(d) and |l7|). For large funnel angles ((3 > 1°), 
the region of increased temperature associated with a 
shock may extend over the intervals between three or 
more packing sites. From Fig. |l8|(a) and(b) we estimate 
that At ~ 2 — 5 cm, hence tt ^0.04 s. 

From plots of T(x, t) it is easy to estimate the decay 
rate of T(x, t) in Eulerian coordinates t t w Xt/U. When 
U » \v m i n \, we find that t t ps t t . Figs, ^(d) and [ij] 
yield the estimate t t % 0.05 s. Generally, it seems that 
t t is always significantly shorter than the characteris- 
tic time separation between neighboring shocks and thus 
fluctuations in the flow velocity created by one shock gen- 
erally do not reach the next shock. 

A more extended region of increased T(x,t) is often 
observed where new shocks are created (see Figs. 0(d) 
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and [l7]). This may be due to the less efficient packing 
of balls in these shocks (which may be related to the 
inefficient energy dissipation and large values of |u m <n| 
discussed in Sec. 0). 

In the above discussion of magnitudes and decay rates 
of T(x, t) it should be recalled that T(x, t) is based on 
measurements of v x and v y with time resolutions ~ 0.02 s 
(see Sec. II C). Granular temperature associated with 
very high collision rates therefore cannot be resolved. 
Due to these high collision rates, the unmeasured part 
of the granular temperature must decay even faster than 
the measured part, and consequently better time resolu- 
tion would merely increase the magnitude of T(x, t) in 
the central shock region and not change the observed de- 
cay rates. 



VIII. AVERAGE PROPERTIES 
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FIG. 19. Spatial dependence of the time averaged relative 
density p(x) showing the packing site periodicity and decay 
near the outlet for different funnel angles: (a) (3 = 0.15°, (b) 
P = 0.3°, (c) P = 0.6°, (d) p = 1.0°, (e) P = 2.0°. (All curves 
are composed of two or three data sets, hence, there are some 
gaps in the data.) 

Moving beyond the flow behavior associated with indi- 
vidual shocks, we will now study the spatial dependence 
of the time averaged values of p(x,t) and v(x,t) which 
we denote as p(x) and v(x) respectively. It is important 
to remember that the measurements only cover a time 
span corresponding to 3-15 passing shock waves. It was 
established earlier that individual shock waves strongly 



influence p(x,t) and v(x,t). Thus, measurements with 
only a few shocks (usually intermittent flows) may ex- 
hibit large fluctuations from the "true" mean values. By 
combining the averages of these two or three measure- 
ments from contiguous sections of the funnel, we obtain 
the statistics in the lowest 110 cm of the funnel in the 
measurements presented below. 

In Fig. [l^ we show p(x) for various intermittent and 
dense flows. In the intermittent flows in Figs. |l^(a) and 
(b), we see that p{x) has a slow growing trend through- 
out the observed region, with weak packing site related 
oscillations. For the dense flows in Figs. |l^(c-e), p{x) is 
nearly unity except near the outlet, and packing site ef- 
fects are more pronounced. The general decrease of p(x) 
near the outlet is most likely caused by the passing of 
fewer shock waves there. 
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FIG. 20. Spatial dependence of the time averaged ball 
velocities —v(x) showing weak packing site periodicity and a 
decreasing overall dependence on x. (All curves are composed 
of two to three data sets, hence, there are gaps in the data.) 



In Fig. |20j the average velocities are shown for differ- 
ent values of /3. For f3 < 0.4° the statistics are not quite 
good enough to ensure that the curves fit together. In 
all curves a general decreasing trend is observed, on top 
of which there are weak packing site related effects. Av- 
erage ball speeds are lowest at the packing sites. Since 
the average flow rate is constant throughout the funnel, 
then \v(x)\ ~ l/w(x)p(x). For large /3, p(x) is essentially 
constant and thus \v{x)\ ~ l/w(x). 
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FIG. 21. Schematic drawings of the (a) flow/speed dia- 
gram v(q) and the (b) fundamental diagram q(p). 



IX. COMPARISON WITH TRAFFIC FLOW 



Three flow types have been identified in traffic flow, 
namely, uncongested flow (steady flow of vehicles with 
speeds ~ 100 km/h and low to moderate densities), queue 
flow (slow flow from 0-20 km/h and near maximum den- 
sity), and queue discharge (accelerating vehicles leaving 
a queue flow situation) [|13[ . Traffic data are usually pre- 
sented as v(q) (the speed/flow relation) or as q(p) (the 
fundamental diagram), where v, q, and p are the velocity, 
flow rate, and density, respectively. As shown schemati- 
cally in Fig ||^, both v(q) (a) and q(p) (b) have regions 
associated with the three flow types mentioned above. 
The fundamental diagram often resembles an inverted 
parabola [l^-[l~^1 , which may contain a discontinuity near 
its maximum as shown in Fig. pT](b) . 

In Fig. ^2] we show grey scale histograms of (q, v) and 
(p, q) for four different values of (3. The average values, 
corresponding to v(q) and q(p), are shown as solid lines. 
(Note that we use the relative density p instead of p.) 
As noted in Sec. Ill C| , at any given time the balls in our 
system are either nearly motionless in a shock wave (cor- 
responding to queue flow) or slowly speeding up (corre- 
sponding to queue discharge). No regimes corresponding 
to steady high speed flow (i.e., uncongested flow) were 
found for (3 > 0.05°. 
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FIG. 22. The (q, v) and (p, q) distributions for four different funnel angles: (a,b) /3 = 0.075° , (c,d) (3 = 0.20° , (e,f ) j3 = 0.50° , 
(g,h) (3 — 1.50°. The solid lines are the averages, indicating the flow/speed diagram v(q) or the fundamental diagram q(p) as 
appropriate. 
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For pipe flow {(3 = 0.075°) Fig. ||(a) shows a small 
queue region and then a large region similar to queue 
discharge behavior. In the fundamental diagram in 
Fig. ||^(b) low densities dominate. The q(p) curve in (b) 
(solid line) has a parabolic shape which peaks at p ~ 0.5. 
The data set from which the distributions were derived 
contained only one shock wave, hence, certain regions 
contain no data since some flow states (e.g., v ~ 80 cm/s) 
simply did not occur during the measurement. 

For intermittent flow (J3 = 0.2°)Fig. p2|(c) displays a 
longer queue flow region than Fig. |22](a). It then changes 
smoothly into a relatively broader queue discharge re- 
gion. The (p, q) distribution in Fig. g2](d) is dominated 
by moderately large densities (p > 0.5) and q(p) peaks 
at a slightly higher value (p ~ 0.7) than pipe flow. 

For the denser flows shown in Figs, ^(e) and (g), it be- 
comes increasingly difficult to distinguish between queue 
flow and queue discharge. The queue discharge regions 
are found at lower values of v and q, which is consis- 
tent with the general trend of falling average velocities 
VIII ) and average flow rates (q ~ /3 -0 4 for 
The corresponding fundamental diagrams 



(see Sec. 
(3 > 0.5° 
in Figs. [2 



2jf) and (h) again have parabolic shapes which 
peak at increasingly higher densities (p > 0.8). The 
shapes of the distributions close to p — 1 are not very 
well-defined. We know from Sees. [TvTj| and [VIl] that 
virtually all density variations are related to packing site 
effects at large (3. Consequently, it should not be ex- 
pected that the distribution for p > 0.9 can be resolved. 
For p < 0.9, the behavior in Fig. |2^(h) is most likely 
caused by the manner in which the dense flow leaves the 
outlet. 

In real traffic, shock waves are generally slower than ve- 
hicle speeds (see |l5"| ) and the densities are low (p < 0.3) 
during uncongested flow. These properties are similar to 
the flow conditions in true pipe flow ((3 — 0°), whereas 
our few measurements (of duration ~ 8 s) do not have 
the time and space resolution necessary to catch the full 
range of the shock dynamics. Thus we cannot display 
ball flow behavior corresponding to uncongested flow. 

Flows for (3 > 0.2° are strongly affected by packing 
sites and thus by the changing number of lanes as illus- 
trated in Fig. EL Several consecutive reductions in the 
number of lanes in traffic are seldom found, but a com- 
parison would be interesting. A more general compari- 
son between the manner in which granular shock waves 
and traffic jams propagate in space and time would also 
be interesting, but unfortunately there are practically no 
available v(x,t) measurements in traffic with the appro- 
priate temporal and spatial resolution. 



X. SUMMARY 

We have presented the results of particle tracking mea- 
surements on a two-dimensional flow of monodisperse 
balls. The circumstances of the creation and propagation 



of shock waves have been studied. In particular, we have 
observed how individual balls behave both in the shock 
region and between shocks, and established that the in- 
teraction of balls between shocks can be disregarded by a 
simple rescaling of gravity. The simultaneous tracking of 
thousands of balls has been used to study the time and 
space variations of quantities such as flow velocity, accel- 
eration and granular temperature. This has been used to 
study the shock wave behavior in very dense flows, the 
processes surrounding the creation of new shocks, and 
the role played by granular temperature. 
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APPENDIX A: BALL POSITION 
DETERMINATION 

Each stored video frame is recalled from the harddisk 
as an array of unsigned 8-bit integers. This array is equiv- 
alent to a grey scale image where black corresponds to 
and white to 255. Each frame is subtracted from the 
image of the empty funnel. Thus, the balls stand out as 
a white patch on a black background. We will call these 
patches ball images. Based on the local distributions of 
grey scale values, a threshold is determined. Everything 
above the threshold is defined as part of a ball. Areas 
of adjoining pixels with values over the threshold belong 
to the same ball. An estimate of the ball center position 
(x c ,y c ) est (based on a weighted average of the pixel val- 
ues) and height (intensity in grey scale) is made. The 
height and estimated positions (x c ,y c ) es t will be used as 
starting values in a fit of a ball's image to a rotationally 
symmetric "bell-shaped" function representing it. 

Next, we need to determine which pixels belong to the 
ball image and should be used in the fit. An area of 
pixels surrounding the estimated ball center (x c ,y c ) es t is 
selected. Let d be the distance from (x c ,y c ) es t for each 
pixel and let r be the ball radius. Pixels with d < 0.9r 
are all selected. Pixels with 0.9r < d < 1.45r are selected 
if the next pixels away from (x Cl y c ) es t are not higher 
than the value of the pixel in question (if they are, it 
indicates that the value of the pixel is influenced by a 
wall or another ball). All other pixels are disregarded. In 
this way, most of the pixels belonging to the ball image 
and some of the surrounding black ones will be selected, 
while pixels belonging to other balls will be left out. 

Typically, between 20 and 35 pixels are selected for 
the position fit of each ball. These pixel values are fitted 
to a rotationally symmetric bell-shaped function with 4 
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fitting parameters: peak height, (x, y) components of the 
center position, and the ball radius. (The convergence of 
the fits is more stable when the radius is allowed to vary.) 
When a radial distribution of pixel values is found, the 
resulting curve matches the function e~' d ' remarkably 
well, where d is the pixel distance from (x c , y c ). In most 
fits, a bell-shaped function with a radial dependence ~ 
e~l d l also gives the best results. In a small number of 
fits (< 1%), however, a stable fit is not achieved with 
this function. In these cases a second fit is made with 
an ordinary Gaussian function e 11 . The Gaussian 
fits give slightly inferior results on average but are much 
more stable in the "difficult" cases. All fits were made 
using the amoeba method At this point a table of 

ball center positions has been determined for each frame 
in the sequence. 

The error of the ball center positions (< 0.15 mm) 
was determined in two ways. In the first, a number (~ 
20) of balls were glued to a rigid, transparent piece of 
plastic. This piece was translated and rotated in the 
funnel and a number of frames were shot. In each frame 
the relative positions of the balls should be constant and 
the noise level could subsequently be determined. In the 
second, individual balls were allowed to roll freely in the 
funnel (without touching the funnel walls or other balls). 
^From the smoothness of enlargements of the resulting 
trajectories, the noise level of the position detection could 
be estimated. 



APPENDIX B: INTERFERENCE WITH PIXEL 
ARRAY OF CAMERA 



In original versions of ball center position histograms 
of the type shown in Fig. [|, there were clear signs of hor- 
izontal stripe patterns where the distance between the 
stripes corresponded to the distance between neighbor- 
ing lines of pixels in the camera. This indicated that 
there was some interference between the pixel periodic- 
ity of the camera and the determined ball center posi- 
tions since there is no reason why some positions relative 
to the pixel array should occur significantly more often 
than others. Subsequent histograms of the fractions of 
the ball center positions in pixel coordinates were made. 
An example of such a histogram is shown in Fig. p3|(a). 
These histograms confirmed the existence of interference 
along both the x axis and the y axis. For simplicity, 
we assume that this interference was the same in all ar- 
eas of the camera and that this systematic remapping is 
"smooth" . With these assumptions it is relatively easy to 
map the positions back to their "true" values and thus 
eliminate the effect of this interference. The resulting 
repositionings of the ball centers are less than half the 
statistical uncertainty of the positions and thus it is of 
limited importance. In Fig. p3|(b) a comparison is made 
between the same ball trajectory before and after the 
remapping. As can be seen, the effect is relatively small. 



The remapping has nevertheless been done for all data 
sets to avoid any "cumulative" effects of this systematic 
error. In Figs. |](a,b) there is still a faint set of horizontal 
lines (with a vertical periodicity of ~ 0.6 mm) but this 
has been significantly reduced by the remapping. 
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FIG. 23. (a) Distribution of fractional part of ball coor- 
dinates measured in pixel coordinates (f3 — 0.5°, D = 10mm, 
6 = 4.1°). (b) Comparison of v y for the same ball before and 
after re-mapping of ball coordinates. The "after re-mapping" 
curve is displaced —2 cm/s. (/3 = 0.15°, D = 10mm, 
= 4.1°). (The ball encounters a shock at t = 0.1 s.) 



APPENDIX C: ASSUMPTION OF TRANSVERSE 
UNIFORMITY 

Throughout this article the flow of balls has been 
treated as essentially one-dimensional. Since the width 
of the funnel (< 10 ball diameters) is relatively small 
compared with the other important length scales in the 
system (e.g., funnel length and shock separation), it is 
reasonable to assume that the flow is uniform across the 
funnel. Between shocks, balls essen tially move indepen- 
dently of each other (see Sec. Ill C ) and uniformity may 
be harder to maintain. In any case, this assumption 
should be checked. It is unfortunately hard to devise 
any meaningful statistics that could confirm the instan- 
taneous transverse equilibrium of a moving shock wave. 
By averaging over x and t on the other hand, we can 
achieve meaningful statistics and subsequently check the 
y/w(x) dependence of various flow properties. (When av- 
eraging over x is done we must renormalize y with w(x) 
to make the averaging meaningful.) 
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FIG. 24. The space and time averaged velocity profile ver- 
sus y/w(x) for different funnel angles (D = 10 mm, 6 — 4.1°). 

In Fig. |i|the y/w dependence of v is shown averaged 
over both x (in 10 cm intervals) and t for various flows. 
Both pipe flow and dense flow show flat velocity profiles 
while intermittent flow ((3 — 0.15° and (3 = 0.3°) exhibits 
a drop near the funnel walls. This is not evidence of 
a shear flow situation. It simply means that between 
shocks (where high ball speeds occur) the regions are 
partially "statistically empty" , namely, when balls leave a 
shock, they uniformly occupy the full width of the funnel. 
The balls closest to the funnel walls will often collide with 
the wall and move towards the center of the funnel. In 
the middle part of the funnel they will collide with several 
other balls and in each collision the transverse energy 
will be reduced and the transverse momentum will be 
averaged out (since balls are coming from both walls). 
As a result the density will be somewhat higher in the 
middle of the funnel compared with regions close to the 
funnel walls. Thus, in the overall statistics, the fastest 
balls mostly contribute in the central part of the funnel. 
In pipe flow (/? = 0°) this effect is not seen because the 
density is too low to eliminate transverse energy in the 
center of the funnel, since each ball basically follows its 
own "zig-zag" trajectory. (Note that the data for (3 = 
0° has been shifted downwards for clarity. The actual 
velocities are ~ 120 cm/s.) In dense flows (/? = 0.6° and 
/3 = 2.0°) the effect is not observed since the density is 
constantly high (see Figs. |l9|(c) and (e)) and there is no 
space for balls close to the funnel walls to move away as 
they accelerate following a shock. 
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